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i.  ijrnionucrior? 


In  1046,  to  estimate  the  location  parameter  0  of  normal  distribu¬ 
tion  N(9,l),  Hosteller  considered  estimators  of  t!io  form 

1c 


0=  l 
i=l 


C-'“(ni), 


(1.1) 


where  . .  <_  are  the  order  statistics  of  a  random  sample 

Xi . y  from  N(6,l),  m  =  [nc^J  +1  with  0  <  -;...<  cy  <  1 ,  Cj's  are 

constants  and  [nrx.]  is  the  largest  integer  less  than  or  equal  to  nno  . 

i-  i  ‘ 

Assuming  equal  weights  c^  =  ...  =  c  =  1/k,  for  each  1;  =  1,2 . he 

found  the  optimal  choice  forio^, . . .  ,a„_)  for  which  the  asymptotic  vari- 

ance  of  0  is  minimized.  In  this  note,  we  generalize  his  work  into  two 
directions.  First,  rather  than  normal  distribution,  an  arbitrary  dis¬ 
tribution  of  a  location  parameter  0  with  pdf  f(N-f)  is  considered,  where 
f  is  assumed  to  be  known.  Second,  instead  of  giving  equal  weights  Tor 
(Cj, . . .  ,c^) ,  the  opt 'ran  weights  n"o  d/dornino''  along  with  the  optimum 


spacing  for  (n^ . n^)  or  ( ,...<%) .  Although  our  argument  is  appli¬ 


cable  to  general  =  1,2 . because  of  analytical  corpl  i<  a*  ion  and 

because  of  satisfactory  achievement,  of  relative  -'sy  intot  ie  of  fie  iency, 
wo  discuss  the  case  of  k  =  2  in  details  and  briefly  treat  the  ease  of 
I;  »  5.  In  section  2,  the  result  js  applied  to  i  be  Cauchy,  distribution 


with  pdf 

f(v-O)  -  l/TTll+Cx-'-)")  ,  -  '  •  (1.2) 

and  the  logistic,  distribution  with  pdr 

■•.--.ro  o 

-  d  -  (1.2) 

as  well  as  normal  distribution.  As  is  veil  Itviwn,  a  omiputat  ional  Jif- 
ficulty  is  involved  jjn  ilerivinq  the  'If  ’s  (mavii  in;  1  i  la  1  ihc  >■  . '  f  t  i  ran  tors) 


f (x-0)  =  c'(X  e)/  (l+e"f>:-ft) 


of  0  for  the  distributions  (1.2)  and  (1.2).  On  the  other  ban.',  it  is 
also  well  known  (e.g.  1.11)  that  for  any  f(x-8)  satisfying  certain  regu¬ 
larity  conditions,  with  1;  =  n  and  suitable  r  =  c.(n,f)  (hnown  as  score 

A 

function),  0  in  (1.1)  becomes  a  DAT!  (best  asymptotically  normal)  rstim- 
ator.  Compared  to  the  HAI-T  estimators,  our  esi  ir.iat.ors  have  a  great  com¬ 
putational  merit  and  yet  are  nearly  efficient.  In  fact,  the  esi inn tors 
in  this  note  use  only  0  or  5  order  statistics  out  of  n  order  statistics 
however  large  n  may  be  and  attain  at  least  Of)’.',  asymptotic  efficiency, 

t 

relative  to  the  CAN  estimators.  More  specifically,  in  the  Cauchy  case 
the  RAT!  {relative  asymptotic  efficiency)  of  the  median  alone  is  Ml'", 

I?ut  the  RAT's  of  our  optinum  estimators  with  k  =  2  and  \  -  5  are  27.' 
and  OCd  respectively.  In  the  logistic  case,  while  the  RAT!  of  ihe  median 
alone  is  ortly  75ft,  the  RAR  of  our  optimum  estimator  wit):  1;  -  2  is  sur¬ 
prisingly  02. 8f';.  It  is  noted  that  the  optimal  weights  r.’s  in  these 
estimators  are  constants  independently  of  n.  lAudher,  ow  estimators 
are  briefly  compared  to  similar  estimators  considered  in  c.ast  wirth 
(1906). 

Wo  remark  that,  associated  with  robust,  estimators,  estimators  of 
the  form  (1.1)  with  k  -  n  have  been  treat  in  groat  many  papers  (e.g. 
LSI). 


2.  Main  result.  Tot  7  bo  the  class  of  coni  iivvi-dy  .'i  T.-vont  ir’de 

,  .  1 

ana  sjfV.t.Kitric.  pdi  y  on  :  with  respect  to  t’  <  T<  ’  law:..  •  wsure,  m.d 
suppose  ICj . :C  are  i.i.d.  wit:-,  pdf  r(*-- n)  •..hen-  f  t  y  i  ;  ’.i,  c-i  and 

"•u'l  '(.)• 

for  any  h,  let  0  <  <  . . .  <  <y  <  1  !-.e  rivd  mbr  an’  lot  f  t'e 

-  quart j le  of  the  population  when  n  =  oj., t 


let  R.,N<  ...  <  7,  .  be  the  enn-es;  viii-’  i  iv  i  .  ,j 

( 1 )  -  -  ( n ) 


"U“ 


"  f 1  r(;<Kx  -  !•(.  .) 


CM) 


’.Vo  assume  that  f(£^)  >  0  (i=l ..  ,k).  'Vrine  tlx*  sample  quantiles  by 

X  ’r  where  n.  -  [nro  1  +  1  (i=l, . . .  ,k).  To  esl.  hinto  the  location 

'ni'  1  1 

parameter  G,  v.t:  consider  the  class  or  linear  mt ir-aior;-.  of  v(,n  of 

i ' 

the  iom, 

*  1 
0  =  OjT.^  j+.**+  a,."^  v  whore  . ;>,)  r.  M.  (•'••’) 

1 

As  is  well  known  (Cl  1  p.  17),  as  p  -►  <N>  joint  ILstriNitinn  of 

t''(nl  )*  ”  ’  ’**(n  onnv('I'^Ki  to  a  k-vnriato  norm]  did  rihul.  ion  'vith 

mean  (0  +  £,,...  ,0  +  f_  )  and  covariance  mlri:;  ~((n.  .)),  where 
-t  j  \  n  ],) 

otj  =  «id-aj)/f(i;i)r(r..)  (i  <_  j). 

Therefore,  for  any  0  in  (2.2),  the  asymptotic  did  ri hut. ion  of  &  i.s 

k 

normal  with  mean  u  «  l  a  (Q  +  £, )  on:';  variance  a  /n,  vl-ore 


1  »i"i (l-ap/rup  <•  l  ■vV'(0-.jVrci 

i«l  1  1  1  1  i  <  j  <’  5  -1  1 


>r(C.) 

I  .1 


1  >  .1 


(2.4) 


Tor  0  to  lx:  consistent  for  0,  t.ho  conditions 

1:  k 

J  a,  «*  1  and  ?  n . C 

ii\  1  jfi  11 


c.r>) 


arc  necessary  ans  sufficient.  let,  C,.  he  t.he  da:  ,  of  e;;t  k  inters  of  the 
form  (2.2)  satisfying  (2.5).  Then  findin;;  n  rnini'mn  asymptotic  variance 
estimator  in  C.  is  equivalent  to  minimi  Tin;;  (2.4)  if  respect,  to  (a  ... 

Ik  1 

a^)  and  (o^ . Oj.)  subject  to  (2.5).  Mo  to  the  retry  of  f,  for  k 


odd,  wo  include  the  met  lion  ^  ond  choose  the  oilier  rjunn1.il  os 

symmetrically  alx>ut  the  median,  i.e. ,  r^  =  l-cy,  =  1-r^  .  or 

equivalently  ^  ^  =  -^_1 , . . . ,  and  the  weights  =  a,. ,  = 

Vi  •  •  •  •  .  Then,  clearly  the  second  condition  of  (2.5)  holds  and 

fUi}  =  f(y>  f(^)  =  . 

Case  k  =  3.  In  this  case,  0  in  (2.2)  is  of  the  fora 

where  •ind  ^ro  the  sa/qile  and  (1-a^)  quantiles  remeotively 
and  X(2)  is  the  median.  From  the  first  condition  of  (2.r>),  "a+’^l  . 

’Writing  f(^)  =  fj  and  f(f>2)  =  f(fl)  =  f o ,  the  asymptotic,  variance  of  0 
in  (2.4)  is  computed  as 


2a,  4a 


(a.a  )  =  a~l— ^  +  -L|  +  al— L  -  +  I— 

f  “  1llo  r  Vo  C  -ir 

t  o  o  o 


(2.7) 


Minimizing  this  with  respect  to  a  yields 

i  2ol  2di  'Iff 

a*  =  a*Caj)  -  r-j-  -  YJ  J/2r- TT  -  j-r  +  "or 
fo  l  o  v  c; 


(2.8) 


O'^a^.aj)  =  a1(l-2a1)/2[2a]f^  -  4n  f  f  +  fj  J  .  '  (2.5) 

To  further  minimize  (2.0)  with  respect  to  (n  <  f(  <-1/2),  let 


0(a,)  =  fx/to  -  f(F  ](a1))/f(0)) 


(2.10) 


where  7  is  the  inverse  of  F  in  (2.7).  ^ert  mjninizinrr  (2.n)  is  equiva¬ 
lent  to  maximizing 


i/r  o^.aj)  -  (n-^r/njd-?^)  +  2,  (a.io) 

which  results  in  the  equation 

n'Cdj)  -  Cd-lojVaa^i-art^l  n(«1)  =  l/d-Po^)  (2.11) 

In  any  specific  problem,  we  first  solve  (2.H)  for  n^(o  <  a  <  1/2), 

and  then  compute  a+  via  (2.8).  Finally  frm  (2.0)  the  minimum  asymptotic 

o  o 

variance  o"  =  o^a^Oj*) ,  a^„)  is  obtained,  where  is  the  solution 

of  (2.11). 

n 

As  a  remark,  in  order  that  o*  1-e  less  t'-nn  1>'e  nsynrdntie  variance 

n 

of  the  median  X.0*,  i.e. ,  l/4f“  ,  it  is  easy  to  see  that 

v*'/  O 

(fi*’?Vi*)"  *  ° 

is  necessary  and  sufficient.  Fence,  unless 

f(F"1(a1*))  =  2fou^,  (2.13) 

0^  =  e(a*,  ,  «1#)  is  asymptotical ly  better  than  the  -nodi an 

Since  (2.13)  is  implied  by 

f (>:)  •  2f(0)7(x)  for  all  x,  (2.1d) 

it  follows  in  particular  that  in  the  case  of  the  'nu'.le  exponential  den¬ 
sity  f(x)  *  cxp(-|x|)/2,  0„  cannot  improve  over  v^n y  Further,  we  remark 

that  one  of  the  optimal  weights  (a*,b*)  may  t.a’o  ivrtat  ive  value. 

Case  k°S.  In  this  case,  0  in  (2.2)  is  of  the  form 

0  -  0(a,b,c,o^,<x,)  -  aX(1)  +  b"(0)  +  +  ’"M)  +  •'1-"(r,v  (n-]r,'> 

•v» lero  ar>r>  are  respectively  ihe  n,  ,  (1-0, 

and  (1-Oj)  quantiles,  and  is  the  median  (»>-  (  l /").  "he 


first  condition  of  (?.f>)  yields  2a  +  2b  +  <■  =  !.  Id  1'^  - 
f^(L)  »  f„  and  f(0)  »  f  .  Then  o"  in  (2.4)  ro-’ucor-,  to 

o/''(a,1t,c,ct^,a^)  «  a 'A  +  h"'-'  +  abr.  +  a1'  +  +  l/4f“ 


where 


2c‘l  4oi  1 

A  «  — rr  -  r— f  +  ~rr  , 

f  “  Vo  f" 

1  o 


2a0  4  a,,  , 

r  53  ~!T  ~  77 

f  n  r 

2  '  o  * 


r  = 


4n  4(v  4*v 

<<  -j  o  '■> 

*■  *  77  ~  r"  r  ~  ?>'  +  Tf 

12  lo  ?'o  f 

o 


p  * 


f,r 

1  o  1 

o 


/ 


(2.17) 


Vinimizinn  (2.10)  with  respect  to  (n,’->)  yiol  ’• 

a#  »  (nr  -  i’^)/('ir  -  V 

\  a  (OAT  -  V>)/(ir  -  i:r) 

and 

o 

o  rp-  +  T','ir'  _  a i 

o“(n*  b*.  f!„  a  ,  c^)  =  -Vr— ^ - — -  +  -V 

!'•  -  4 AC’  1!"' 


(7. 10) 


(2.10) 


For  ready  reference,  we  record  below 


/I 

i  r 

4o0-l 

'f2 

-  7“ 

o' 

lbro 

1‘Vtj  n  i 


ba. 


10a,  re 


-t-  1  ,alr  _  _ 

S  THr  Vo  riro 
In  •<  o 


1 

v-> 


1/1  IV 


i  o  r 


Vo  V 


and 


C'“  +  m;  -  = 


^  i_::  j: 

,riro’fn  ii: 


v,d  r. 
^  r 


\2  [Ca, 


r  r  . 

2  o  f 


.Vo  '  r2  “  r 

o  1 


,r’l  /I  1  \  p 


(2.21) 


l/i  i  \  j  _1  .  I  l_I  JI 

7'7n  rJ  lV«  r.  fr«  rj- 


Hxact  algebraic.  minimization  of  !'A,  f'+,  .  b,)  '''ll1'  respect 

to  and  seems  difficult  althcniqh  in  specific  nrnMms,  the  solution 

can  1>e  obtained  numerically.  It  is  remarked  mb  ovon  an  approximate 

» 

solution  of  (cl,  «n)  coupled  with  the  optimal  woipht.s  (a*,  b*)  in  (2.12) 
may  lie  useful  in  Getting  0  letter  than  the  nodian  and  closer  to  Hiernoff, 
Ga.stwirth  and  .Johns’  ’’.A?’  os  t  inn  tor  (son  n  '). 

G.  Comparison  of  ost  inators.  "lie  table  l  hole-  .•nr.na.rir.es  feature*-- 
of  our  optimum  estimators  for  k  =  2  for  normal ,  Cauchy  an  1  logistic. 

distributions. 


Normal 


i  Qiuchy 


I  opiatic  I  0.2500  O.G 


a* 

V>+«l-2u*  j 

-i 

'•s 

V 

0.3054 

0 . 3 SO 2  I 

1  .non 

1 

nrf  j 

(0.3) 

(0.4)  | 

(  ^  .  eo;jf ,) 

(.7 ?>"■  ) 

-0.0207 

• 

1.0114  j 

o 

so.fr/, 

(0.3) 

(0.4)  ! 

(  •  >  _  M  lp  ) 

0.3 

n.4  j 

p 

op  pe 
•  * 

(0.3) 

i 

(0.4)  j 

(3.:nr,7) 

!  (rrr. 

Here  the  numbers  in  the  parentheses  are  those  rfriiir’  f ron  the  view¬ 
point  of  robustness  by  Gost.wirth  (1000) ,  and  1/T  is  f'e  asy-atiit  jc 
minimum  variance  of  tlie  NAN  estimator  where  T  denotes  l  ho  risiv'p  infor¬ 
mation.  In  the  normal  case,  while  the  r''T"  of  ‘Vv.-tol  s  «-;•:»  i-ml'T 
with  k  m  G,  =  0.1320  and  equal  wpirhts  is  PfT,  the  ”.v  nf  our  estim¬ 
ator  is  DOC,.  Corpare-d  to  Onst.’ri rth ' s  robust  estimator,  "-hi  eh  is  V'o 


same  fom  ru^  our  ost inators ,  naturally  our  v.,omtors  perform  l  etter 
since  in  our  case,  knowing  the  form  of  f,  the  optimal  woi'du.s  with 
optimal  spacing  arc  obtained.  This  might  suggest.  that,  when  n  is 
large,  it  will  be  1  letter  to  try  to  know  the  fom  of  f  rather  than  to 
apply  "robust”  procedures.  However ,  for  example,  test. in;:  the  normal¬ 
ity  in  an  efficient  way  is  not  an  easy  task.  Tinco  for  k  =  2 ,  the 
RAH  of  Cauchy  case  is  still  less  th;m  ft V-,  w  computed  the  ease  of 
k  «=  5  numerically,  hie  resulting  optimum  values  are  as  follows 

V 

a*  =  -0.0720,  h*  =  0.0135,  c*  =  0.0312 

ct^+  =  0.10,  an±  =  0.10  and  o'*  =  2.22  . 

In  this  case,  the  PAR  of  0*  is  00.]','  .  It  is  noted  i  bat  wh.en  U  =  3 
and  k  =  5,  some  of  the  optimal  weights  are  negative,  which  is  consis¬ 
tent  wit’-,  the  result  by  C-.emoff,  ra:?t-irtk  a-/  .To*  r.  ■' ' . 
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